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Abstract 

In the context of the dynamical evolution in a non-stationary thermal bath, we construct a 
family of fluctuation relations for the entropy production that are not verified by the work 
performed on the system. We exhibit fluctuation relations which are global versions either of the 
generalized Fluctuation-Dissipation Theorem around a non-equilibrium diffusion or of the usual 
Fluctuation-Dissipation Theorem for energy resulting from a pulse of temperature. 

1 Introduction 

One important recent progress in non-equilibrium statistical physics is the discovery of var- 
ious fluctuation relations which can be viewed as non-perturbative extensions of the usual 
Fluctuation-Dissipation Theorem (FDT) |29ll36| . Such relations pertain either to non-stationary 
transient situations |141 125] or to stationary regimes [16]. In particular, a family of fluctuation 
relations holds for the distribution of work performed on a system |25l [TOl H5l [5] which evolves 
in an equilibrium bath. This developments had an important impact on the physics of nanosys- 
tems and biomolecules j43]. Here, we examine the question of the extension of such relations 
to an evolution in a non-equilibrium medium. More precisely, we consider systems placed in 
the thermal bath that is non-stationary, a situation which can be realized experimentally by 
adding or extracting heat, or by modulating the pressure. Let us remark that non-stationarity 
of the temperature is linked via Stokes law to non-stationarity of friction and that particles with 
time-dependent radii have non-constant mass and friction. This last situation is an important 
problem in astrophysics for the formation of planets through dust aggregation (11 [l]. There is 
another situation where the non-stationarity of friction can be realized with particles diffusing 
in ferroelectric fluids when external magnetic fields are controlling the intrinsic viscosity [22] . 

A famous example of a system which evolves in a non-stationary bath is the Temperature 
Ratchet Model [H] of a Brownian motor. We can also consider the evolution of a system initially 
in equilibrium at high temperature (i.e. in the state with the probability density pi = exp(— |^)) 
put in contact with a non stationary bath reaching a low temperature Tf. The function T{t) is 
then called the cooling schedule. A particularly interesting case is the instantaneous quench of 
the initial system in a thermostat at temperature Tf which correspond to the cooling schedule 
with instantaneous initial change of temperature. The age (or waiting time) of the system is 
then the time elapsed since the quench. Finally, at the more formal level, such non-stationary 
evolution rules appear naturally from a stationary dynamics after system size or coarse grained 
expansion [19]. Diffusion properties in a non-stationary medium have been recently studied in 
the case of unidimensional overdamped dynamics driven by multiplicative non stationary noise 

[mm]. 

The present paper consists of five sections. Sect. 2 sets the stage and notations for the model 
of non-equilibrium and non-linear Langevin dynamics that we consider. In Sect. 3, we recall the 



main result of (5] and we list different time inversions that permit to obtain fluctuation relations 
for entropy production. We observe that in this context the work performed on the system no 
longer obeys the fluctuation relations but that there still exists a functional that upon averaging 
gives the free energy change in a non-equilibrium protocol. In Sect. 4. we explicitly construct 
fluctuation-dissipation relations that result from the Taylor expansion of the last fluctuation 
relation. In particular, we recover the generalized FDT f6] around a non-equilibrium diffusion 
and the usual FDT for the energy which results from a pulse of temperature |42j. We discuss also 
the physical meaning of the effective temperature of an harmonic oscillator in a non-stationary 
bath. Finally, Sect. 5 presents our conclusions. 

Acknowlegements. The author thanks Krzysztof Gawedzki, Frangois Delduc and Gregory 
Falkovich for discussions and acknowledges support of the Koshland center for Basic Research. 

2 Thermally Driven Non-Linear Langevin Dynamics 

The general dynamics that we consider is described by a non-autonomous diffusion process that 
we call Langevin dynamics (or non-linear Brownian motion) in a d-dimensional (phase-)space E. 
The process obeys the stochastic differential equation (SDE) 

= -Tnx)d,H,{x) + nY{x)d,H,{x) + G\{x) + w\{x) + r;j(x) (1) 

which should be interpreted in the Stratonovich convention. Here Ht{x) is the Hamiltonian of 
the system (the time indexation corresponds to an explicit time dependence), Vtix) is a family 
of non-negative matrices, Ii-t{x) is a family of antisymmetric matrices, Gt{x) is an external force 
(or a shear) and wt{x) is an additional corrective term which will be explicited in eq. ([s} below 
and which comes from the a;-dependence of and lif Finally r]t{x) is a white Gaussian vector 
field with mean zero and covariance 

{rilixHiv)) = 2Dnx,y)5{t-s) with A(x,x) = dt{x) = ]^Tt{x). (2) 

Pt 

The last equality, which contains the temperature of the bath Tt = is called the Einstein 
relation. In the case when the bath is in equilibrium, this equation reflects the fact that the 
friction and the noise are two dual effects of the interaction with the bath. In this non stationary 
setup, the Einstein relation is valid under the assumption of local thermal equilibrium, provided 
the thermalisation of the immediate environment of the Brownian particle is much faster than 
the characteristic variation of temperature. Let us underline that in this context the term "non- 
linear" concerns the non-homogeneous properties of Ft (a;) and \it{x) (which may characterize non 
homogeneous properties of the bath) and not the fact that equation ([ij is non-linear. Such non- 
linear properties appear naturally in many situations: non-ideal plasmas and gases |26ll27] . ultra- 
cold clusters of atoms or molecules cooled by interaction with laser radiation, active Brownian 
particles [13]. Finally, the additional corrective term in ^ is given by the expression 

w\{x) = dy,Dnx,y)\^=y ~ ^a,nj^(a;). (3) 

Pt 

The presence of this term assures that in the case with stationary Hamiltonian and temperature 
(i.e. Ht = H, f3t = (3) and without external force (i.e G = 0), the Gibbs density exp{—PH) is an 
invariant density]^ The presence of this term Wt{x) in equation ([T| can appear as a makeshift 
arrangement, but it was extensively studied in the literature of non linear Brownian motion 
|26l 127) . Note that Wt vanishes in the case of linear Brownian motion where Dt{x,y) = Dt and 
Ht{x) — Hf. We call the deterministic part of the second member of eq. ([TJ: 

ut{x) = -Ty{x)djHt{x) + Uy{x)d,Ht{x) + Glix) + wl{x) (4) 

^In fact, the Gibbs density is then an equilibrium density: an invariant density with vanishing modified probability 
current [6]. 
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the drift term. An elementary case of non-Hnear Brownian motion is the Landau-Lifshitz-Bloch 
dynamics of a Brownian spin jl^ in an effective magnetic field B^^^ {x) = —VHt (which can 
incorporate interaction with other spins). It follows the dynamics 

x^-xxVHt + Xtx x{xx VHt) + Gt{x) + x x Ct with {QCi) = -^5''5{t - s). (5) 

Pt 

The first term on the right hand side is the precession term, the second one is the damping term 
and the third one an external torque. The noise (and the damping term) accounts for the effect of 
the interaction with the microscopic degrees of freedom (phonons, conducting electrons, nuclear 
spins...). This dynamics is a particular case of ([l| with T'/{x) = Xt ((5*^ x^ — x^x^ ) and n^-' (x) = 
^ijk^k Tffii]^ ^tjk |.jjg totally antisymmetric tensor. In this example one sees the need for the 
term Hi'' {x)djHt{x) in eq. all corresponding to the Hamiltonian vector field which also permits 
to describe systems of non-overdamped Brownian particles in the phase-space with coordinates 



), non-stationary mass, in an external potential Vt, subjected to a non-conservative force ft 



and in a non-stationary bath giving rise to the noise and the non-homogeneous, non-linear but 
isotropic drag 74(5, p). The Stratonovich SDE which governs this model is then 

which is again a sub-case of ( 1 1 with Ft = ( ^ ^ ) , 11 = ( ) , Ht — -f Vt {q) and 

Gt — { f^(^q) )■ "^^^ Kramers case corresponds to the Stokes law of friction jt{q,p) — lt{q)- 

Another example of friction with 'jt{q,p) = lip^ —Pa) appears in the Rayleigh-Helmholtz theory 
of sound |40] . 

We start by collecting the elementary properties of diffusion processes that we shall need 0^ . 
The Markovian generator Lt of the process xt satisfying the SDE ([ij is defined by the relation 

dt{f{xt)) = {{Ltf){xt)) , where Lt = UfV + V- ^V, (7) 

where the modified drift ut is defined in term of the drift ([i} : 

^\i.x) = u\{x) - dy,DY{x,y)\^^y. (8) 

The time evolution of the instantaneous probability density function of the process pt{x) — 
{5 {xt — x)) is governed by the formal adjoint l\ of the generator Lt 

dtPt = L\pt (9) 

which can be rewritten as a continuity equation (resp. an hydrodynamic advection equation) by 
defining the probability current jt (resp. the mean local velocity Vt^ 

dtPt = -V • Jt = -V • {ptvt) where Jl = (uj - l3t^Tl'V,)pt and vl=^-^. (10) 

Pt 

As was explained in [6], it is convenient to use the freedom to add a divergenceless term in 
the definition of the probability current to obtain the modified current and the modified local 
velocity 

j,'^jl + /3r'V,(ffV0 and vl^^-^ (11) 

Pt 



which verify also the continuity equation ( |10| but vanish in the case of stationary Hamiltonian 
and temperature (i.e. Ht — H , f3t — P) for vanishing external force Gt — 0. 



3 



3 Fluctuation relations and time inversion 



In [5], various fluctuation relations were discussed for arbitrary diffusion processes. We recall 
here the main result in the context of systems with dynamics of type (jlj, see also Sect. 3 of 
[3- With the use of combined Girsanov and Feynman-Kac formulae, one obtains the detailed 
fluctuation relation (DFR) 

fi,{d^) PAx;dy,dW) e-"^ = fil{dy*) {y* ; dx* , d{~W)) , (12) 

where 

1. fJ-oidx) — g„{x) dx is the initial distribution of the original forward process ([ij, 

2. iJ-lidx) = q'^^{x) dx is the initial distribution of the backward process obtained from the 
forward process by applying a time inversion (see below), 

3. Py(a;; dy, dVF) is the joint probability distribution of the time T position x^ of the forward 
process starting at time zero at x and of afunctional (linked to the entropy production) 
of the same process on the interval [0,r] (described later), 

4. (x; dy, dW) is the similar joint probability distribution for the backward process. 
The time inversion acts on time and space by an involution 

{t,x) ^ {t* = T-t,x*). (13) 

Such an involution induces the action x t—^ x on trajectories by the formula Xt — x^_t and, 
further, the action on functional of trajectories F 1-^ by setting F[x] — F[x]. To recover 
various fluctuation relations discussed in the literature |28l ISTl [TQ| [25| HgI |8], one divides the 
drift part Q of ([iJ into two parts, u — u+ + U-, with m+ transforming as a vector fleld under 



the space-time involution ( 13 1 and u_ as a pseudo- vector fleld: 



u"j,_t^^{x*) = ±{dkx*'){x) ul^{x) , u'' = m; + w1 . (14) 

The random fleld rit{x) may be transformed with either rule. By deflnition, the backward process 
satisfles then the Stratonovich SDE 

X = ulix) + ril{x) (15) 

and, in general, differs from the naive time inversion xt of the forward process. The functional 
is given by the expression 

T 

=ln£-„(xj -\n{Ql{x*T)(j{xT)) + Jtdt (16) 







where a{x) = det 



is the Jacobian of the spatial involution. The intensive functional Jt 
has the interpretation of the rate of entropy production in the environment and is given by the 
expression 

Jt = PtUt,+ {xt) ■T~^{xt){xt-ut-{xt)) - {V ■ut^-){xt) , (17) 



The time integral in eq. (16 1 is taken in the Stratonovich sense. When ^^{dx*) — iir^{dx) then 
the boundary contribution [In Qgix^) — ln( gQ{x^)a{xT))] to gives the change in the instan- 
taneous entropy of the process. In this case, the functional W^, becomes equal to the overall 
entropy production. Moreover, with the interpretation of Jt as the entropy production in the 
environment, the First Principle gives us that the work 7^ performed on the system can be 
expressed in term of Jt: 

T 

= HrixT) - Hoixo) + J^dt. (18) 
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We can underline that in this setup, and contrary to the case of a the stationary bath jS], the 
work 7^ cannot be identified with the functional for an appropriate choice of initial densities 
of the forward and backward processes. This means that the work does not verify the DFR. 



The DFR ( 12 1 holds even if the measures /i^ and /i^ are not normalized, or even not 
normalizable. When they are normalized, let us denote by ^ — ) and by ( — the expectations 
of functional of, respectively, the forward and the backward process on the time interval [0,r], 
with initial distributions /i„ and /i^ . One of the immediate consequences of the DFR equation 
(121 is the (generalized) Jarzynski equality [25] 

(e-^r) = 1 (19) 



obtained by the integration of the both sides of eq. (12 I. It implies the inequality (W^) > 
that has the form of the Second Law of Thermodynamics stating the positivity of the average 
entropy production. With a little more work [5], the DFR (12 1 may be cast into a form of the 
(generalized) Crooks relation jlO] : 

(FeM-Wr)) = {Fy. (20) 



We will now restrict ourselves to the class of time inversions ( 14 1 such that there exists a 
non-stationary density ft such that 

u+.t^ P^^TtV In ft and then u-.t ^ ~TtV{Ht + [3^^ In ft) + IltVHt -^V-Uf + Gf (21) 

Pt 

After a straightforward calculation, the rate of entropy production in the environment may be 
expressed as 

Jt = Xf (Vln/t)(xt) + {{ft)-^Llft){xt). (22) 

With the choice fio{dx) — fo{x)dx et iJ.Q{dx*) = fT{x)dx the functional Wt takes then the simple 
form 

Wt = {{ft)-' Lift ~ dt Hft)) (xt) dt. (23) 

We shall see that the fiuctuation relations associated to this peculiar family of inversions are the 
natural generalizations of the FDT. We shall now describe particular cases in this family of time 
inversions. 



3.1 Complete reversal [5] 



As the function ft in ( |2l| we take the instantaneous density function (i.e. ft = pt) of the forward 
process ([ij distributed with initial condition /q. Here, the functional (231 trivially vanishes 
Wt = and the DFR (12 1 takes the form of the generalized detailed balance 



t^oidx) P^{x;dy) = n^{dy) P^{y*;dx*). 



(24) 



One may show that pl{x) = pf{x*) is the instantaneous density of the backward process and 
that the corresponding probability current satisfies the relation 



Jt^'ix) = -{dkx*''){x)j''^,{x*). 

This inversion is employed in many articles in probability theory 
to the vanishing overall entropy production. 



(25) 

. It corresponds 



3.2 Current reversal P [5] 

Another useful choice of time inversion, called the current reversal is based on the choice ft — irt 
where ttj satisfies Lint — —^-jt — 0. In the case where Gt — 0, we have ttj — exp{—Pt{Ht — Ft)) 
with Ft the free energy (i.e. exp(— /3(Ft) = / exp{—f3tHt) )■ One can show |5] that 7rJ'(x) = 
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TTf (x*) is the density for the backward process which correspond to the conserved current V-j[ 
and that eq. (251 still holds. The functional (23) takes now the form 



{dtlinrt){xt) dt, 



(26) 



where the index "ex" stands for "excess" |45ll5]. For the backward process, the functional W':^'^ is 



given by the same expression with tt^ replaced by 7r[. The Jarzynski equality (19 1 for this case 
was first proven in one dimension in |2l] and in the general case in |20l l32l |5]. 



3.3 Canonical inversion [5] 

A natural choice for the system ([ij is to take ft to be the Gibbs density exp{— Pt{Ht — Ft)) where 
Ft is the free energy (i.e. exp{—f3tFt) — J exp{—(3tHt) ) if the Gibbs density is normalizable and 
zero otherwise. This corresponds to the choice w+ t — —TtVHt in (21 1. The functional ( [23) ) 
becomes 



1 

-{PtFt - PoFo) + J [dt iPtHt) + PtGt -VHt-V- Gt] {xt) dt 



(27) 



where the index "ci" means "canonical inversion". The generalized Jarzynski equality (191 can 
be rewritten in the form: 



^exp(- j [dt{l3tHt) + l3r 



GfVHt-V- Gt] (xt) dt) ) - exp [-(/^t^^t - PoFo)] (28) 



which permits to extract the difference of free energy out of a non-equilibrium experiment in a 
non-stationary bath (but the connexion with the work performed is lost). For example, for the 
Brownian particle (|6|, this functional takes the form: 



W^' = ~iPTFT-(3oFo) + 



iQt))+f3t{' 



^ + idtVt){qt) + ift-VVt)iqt)) dt (29) 



which, as compared to the functional which appears in the usual Jarzynski equality |25l [28] . 
contains new terms proportional to the variation of temperature dt(3t and of mass dtmt- From 



the Jarzynski equality (191, one can deduce the inequality which constrains the evolution in a 
non-stationary bath 

(W^") > . 

In the case with a stationary Hamiltonian (i.e. Hf = H ) and without external force (i.e. Gt = 0) 
this constraint reads 

{dtf3t){H{xt))dt > {(3TFT-f3,Fi). 

Let us consider now the overdamped particle in a 3-dimensional time-dependent harmonic 
potential. Such example admits an analytical computation of the distribution of the func- 
tional W^^. We consider a non stationary bath {jt,Pt) and the harmonic potential Ut{x) = 
^{x — at)^ where kt is the stiffness coefficient and at is the instantaneous center of the potential. 
The particle is initially distributed with the Gibbs density Pq{x) = exp{— (3(){Uo{x) — Fq)) = 

3 

(It) ^ ^-'^P ~ '^o)^) • Further, we will restrict our study to the particular case, not nec- 

essarily physical, where the temperature of the bath and the stiffness coefficient are such that 
their product is stationary: ktf3t — C. This setup generalizes the unidimensional stationary case 
{kt = k , fit = P ,lt = 1 and at = ut) considered in j50| 134 j . The system satisfies the linear SDE 



h 

It 



{x — at) + r]t with 



{vlvi)-j-S{t~s)S^- 
Ptlt 



(30) 



6 



The functional W^* — W^^ takes here the form 



W^' ^-C I at- {xt - at) . 
Jo 

The distribution of for a process with the initial Gaussian density po is Gaussian due to the 



linearity of ( 30 1 . A straightforward calculation gives the mean 

i-T rt / rt 



{W^')^c[ dt [ dfd,exp(- [ —dv)ds (31) 

^0 "'0 \ Js lu J 

and the variance of this Gaussian (we assume that the integrals exist): 

V^T = ((W^T - (W^t))^) = 2C^ dt dfdsexp(^~ '^du^ds. (32) 



The distribution of W^^^ is then 



ivj ) '''' 

and an elementary calculus shows that the Jarzynski equality ([l9|is equivalent to the fact that 



Vt = 2 (yV^^ which is evident from the comparison of (31 1 andT32 



3.4 New inversion 

We choose for the function jt the mean instantaneous density p\ of another Langevin dynamics 
([ij which possesses the same parameters Ft, lit, Gt but with another non autonomous Hamil- 
tonian Ji[ and another bath temperature We note Lj (resp. L\ ) the Markovian generators 
of the process with the Hamiltonian Ht and the bath temperature fit (resp. }1[ and /SJ. The 
functional ( [23) ) takes now the form 

Wt = f [(p;)-^ {Lt - L',f p',] (xt) dt. (34) 
Jo 

This new inversion will permit to recover new generalizations of the FDT around non-stationary 
non-equilibrium diffusions, see also |6]|71|3]. 



4 Generalizations of the Fluctuation-Dissipation Theorem 



As noted in |17ll3T] . the fluctuation relations may be viewed as extensions to the non-perturbative 
regime of the Green-Kubo and Onsager relations for the non-equilibrium transport coefficients 
valid within the linear response description of the vicinity of the equilibrium. Ref. jS] contains a 
detailed argument showing that if in a stationary bath one perturbs an equilibrium system by 
introducing a weakly time dependent Hamiltonian Ht{x) — H{x) — ga,tO"'{x) then the Jarzynski 



equality associated to (26 1 or (27 1 gives in the second order of the Taylor expansion in g the 



usual FDT. Ref. [J| showed that similar correspondence still holds around non-equilibrium steady 
states for a stationary dynamics with an external force (i.e. G ^ 0). In this case, it is the Crooks 
relation ( [20| associated to the functional (26 1 which gives the modified Fluctuation-Dissipation 
Theorem [xl [21] after the first order Taylor expansion in g. The second order expansion of 



Jarzynski equality (19 1 associated to the functional (26 1 gives in such a situation only a special 



case of this theorem. We shall now investigate which type of fluctuation-dissipation identities 
may be deduced by Taylor expanding the fluctuation relation corresponding to the time inversion 
of Sect 3.4. 
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4.1 FDT around non-stationary difFusions 

We consider the system with the Hamiltonian Ht{x) — Hl{x)—ga,tO°'{x). Following Sect. 3.4, 
we choose ft as the mean instantaneous density p° of the unperturbed system with 5 = 0. The 
functional (34 1 becomes 



Wt 



r 9a.s {{plr^Mfpl)^ ds with m: = (r.vo'^ - n.vo'^) • v, 



(35) 



where the subscript "s" on ((p°)~^M^^^p^) signals that the latter function should be taken at the 
point Xs, Let us now write a particular case of Crooks relation (20 1, where the average is in the 
system ([ij with the Hamiltonian Ht, associated to a single time functional F[x\ = 0°'{xt) = 01 
{0<t<T): 



(36) 



We shall denote by ( . )g the average of the process with the dynamics driven by and by 
■yO and ti°, respectively, its Markovian generator, its mean local velocity and its modified mean 
local velocity. The first order Taylor expansion 



exp(— Wt) — 1 



(37) 



in ( 36 1 gives the relation 

S 



56. 



9=0 



m ds 



(O? {p°,-Hl'^p^)Jds + Oig') = {Or-tY. (38) 



The right hand side has a functional dependence only on {gu,u > t}, so if we apply 
for < s < t to the last identity, we obtain: 

6 



3=0 



(39) 



A short calculation gives: 



psLlOl + Knd,0')d. In(p2) + {dj\\^)d,0' 



(40) 



= /3.(2«o.V-LO 
= 13s (2«~o . V - L°) O'' - Psv'i ■ WO". 
Moreover, we have the sum rule (for s < t) 



O" - psv^, ■ X/O" + n:^(a,0'')9, ln(p,) + (5,0^^)5,0" 



Of 



With (39 1, (41 1 and(40l, we obtain the Modified Fluctuation-Dissipation Theorem 

S 



1 



9=0 



(on + ((«s-vo^).o,«)g. 



(41) 



(42) 



This is a generalization of the FDT around a non-equilibrium diffusion process in a stationary 
bath {/3t — (3) of refs. [T] |6] [S] [11] and of FDT around non-stationary Langevin equation of 
ref. [52]. It may be also proven as in [6] using the fact that the diffusion process becomes an 
equilibrium one in the Lagrangian frame of its modified mean local velocity w'^, verifying in that 
frame the usual FDT. The transformation of the latter back to the Eulerian (i.e. laboratory) 
frame leads to (42 1. Let us remark that here, similarly as in the stationary case discussed in [7], 
the Jarzynski equality (19 1 for the functional (34 1 leads upon the second order expansion in g 
to a particular case of the MFDT where the observable is replaced by A°' = {pt)~^ M^^ p^ 
which is a (time-dependent) functional of O". 
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The violation of the usual FDT can be parametrized by using (42 1 via the introduction of 
the so-called effective temperature |12l |9] defined by 



((i;0-VO"),O?)^ 



s 



3=0 



(Of) 



(43) 



We shall consider now the case where this effective temperature may be computed analytically 
in order to verify its physical consistency. 



Unidimensional harmonic oscillator in a non-stationary bath. The SDE which 
governs this systems is 



2T, 



1 



x+T], with (rjt-qs) ^ 5{t - s) 



1 



(44) 



We take the cooling schedule Tt such that the bath passes from an initial temperature (Tq — Ti) to 
a lower final temperature {Tr = Tf < Ti) during a time r. The system is initially in equilibrium 
with the bath and its initial density is po{x) = exp{—^x'^)/Zi. We consider two particular 
examples of cooling schedules: 

• Instantaneous quench 



Tf 



if t = 
if t > 



Linear decrease of temperature: 



Tt 



_ r T, + ^{Tf~Ti) ii t<T 



Tf 



if t> 



(45) 



(46) 



Due to the linearity of eq. (|44||, one can compute explicitly the response of the position to an 

V^{x) = V{x) — gtx with go 



external perturbation V{x) 



2"^ 



0: 



1 



g=0 



xt) = - cxp (t - s) 
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k 



1 



(47) 



It has a stationary form and is independent on the cooling schedule. In a similar way, we obtain 
also an explicit expression for the dynamical 2-time correlation function of the position in the 
unperturbed system in the two cooling schedules. For s < i, we obtain for the instantaneous 
quench: 



{XsXt)o 



T, - T 



exp 



^(s + t) ) + ^exp 



1 



and for the linear decrease of temperature schedule: 

{XsXt)^ - 



+ - ^)) exp (-^(t - s)) + exp(-^(. + t)) if 

(-^(* - •^)) - (exp(f r) - 1) exp (-^{s + t)) 



^cxpi 



(48) 



s < r, s <t 
\iT < s <t. 



(49) 

We see in these two formulae that the characteristic time of convergence toward the Gibbs density 
exp(— 2^^^^)/^/ fc foi' the instantaneous quench and t + ^ for the linear decrease schedule. 
Note the relation 



{xsXt)Q = exp ^-^(i- s)^ = (2:^)0 7 



6 

Sgs 



{Xt) 



(50) 



holding for both cooling schedules. It shows that at very large times (i.e. t > s ^ r, ^ ) the 
correlation functions (48 1 and (49 1 take a stationary form depending on t — s. The instantaneous 
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mean density of the process is Gaussian at all times. It follows that the mean local velocity has 
the form 

v°{x) = (-- + ^^)x (51) 



and the corrective term in the FDT (42 I is 



7 



{XgXl) Q 



(52) 



Using the relations (501 and (52 I, the FDT (42 I may be rewritten in this case as the identity 



/ k \ S 
{ds + -) {xsXt)^ = ZTs— (xt) 
V ^/ dgs 0=0 



(53) 



which is easy to check directly. 

The effective temperatures for the instantaneous quench Tq^ and for the linear decrease of 
temperature schedule T^{/ are: 



T^"{s,t,x) = Tf + {Tf 



Ti) exp(- — s) 
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and 



T, + i(r/-T,)(5+^(l-exp(-f5) 



T/ + (T/-r,)exp(-fs) 



cxp(f r)- 



if < s < t 



if < s < r, s <t 
if T < s<t. 



(54) 



(55) 



The two effective temperatures have an expected behavior for large time s converging toward 
Tf. However, we may see in this system the problems with the physical interpretation of the 
effective temperature |12l l9]. For example, HuIs^oTq'^ = 2Tf— Ti ^ Ti and this expression 
can be negative ii Tf < The possibility to find negative effective temperature has been 
observed also in [7] and, for the kinetically constrained model, in |37) . Moreover, the effective 
temperature grows toward its limit Tf, which does not correspond to the physical intuition for 
the temperature of a cooled system. The investigation of the linear-decrease cooling schedule is 
instructive for the understanding of these two problems. For this schedule, lims^o ^li/ — Ti 
and the effective temperature decreases from this value until time t when it reaches T^/J (t) = 



Tf + {Tf-Ti) 



< Tf. So, the problem with the initial time limit of Tg was due to 



the instantaneous modeling of the quench. On the other hand, the second part of the evolution 
for the linear cooling schedule (i.e. for s > t) begins with an effective temperature below Tf that 
may be even negative. The last features are not really physically satisfying but the first one 
explains why the effective temperature Tq^^ converges toward Tf by growing, the fact which is 
confirmed by (551. We represent below in Figure 1 typical joint evolution for the linear cooling 
schedule of the bath temperature T (the crosses) and the effective temperature T^^/ (the solid 
lines) in the case where T, = 300K, Tf = 200K and ^ = Is'^ for r = 0.1s (red) r = Is(blue) 
and r — 10s (black). The singularity of the limit r ^ is evident on this graph. 



4,2 Response of a diffusion to a pulse of bath temperature 

Let us consider the system whose dynamics is governed by eq. ([ij with the variable bath tem- 
perature /3j~^ = {1 + gt)PQ^ . We choose as the function /( the mean instantaneous density of 
the similar system with the constant bath temperature Pq^- The functional (34 1 becomes 



gs{{pT'M^,p',)js with Af, = (/3o)-' ((- 



v,r^)v, + r^ v,v,) 



(56) 
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12 



Time s 



Figure 1: Red solid: r^{/ for r 
Blue crosses: T for r = 



0.1s. Red crosses: T for r 

10s. Black crosses: T for r 



Is. Black solid: T^^/ for r 



0.1s. Blue solid: T^^/ for r 



10s. 



With the same reasoning as in Sect. 4.1, we find the link between the response to a pulse of 
temperature at time s and the dynamical correlation function in the system with stationary 
inverse temperature (3^. 



{0,) = {{{p^)-Hl^pl)^0,{x^))^. 



(57) 



Here, there does not seem to exist a simplification of this relation in the spirit of (42 I and 



we cannot say more except for the case when the system with the bath temperature /3o is an 
equilibrium one. 



Temperature pulse around equilibrium. In the case where the system with the bath 
temperature /3o is in equilibrium (i.e. without external force = and with a stationary 



Hamiltonian Ht — H and the Gibbsian instantaneous density), the functional (27 1 takes the 
form: 



-{PtFt - l3oFo) + / PtH{xt)dt 



(58) 



We want to prove that the Taylor expansion in the second order of the Jarzynski equality p9| 
associated to this functional gives the usual Fluctuation-Dissipation Theorem for the energy |42) 



ds {HsHt)„ - ITT— 

Po ogs 



(59) 
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The equality (19 1 takes now the form 



exp 



(- / ptH{xt)dt) 



J exp {—^tH{x)) dx 
J exp {—f3oH{x)) dx 



(60) 



We develop the left member in second order in gt or ht = gt — gt assuming that gt vanishes for 
t < 0: 

^exp ( - ^ $tH{xt)dt^^^ ^exp (/3o ^ htH{xt) d?j + 0{h^)^ 



2 /-T 



+ / htH{xt) dt - 



PI 
2 7o 



dt / hthsH{xt)H{xs)ds + 0{h*) 



= l + /3o/ ht{H{xt))odt + Po dt hth. 



JO 



3=0 



2 /-T 



+ Y / / htK{H{xt)H{x,))^ds + 0{h') 



"'0 

T 



l + /3o / ht{H{xt))odt + po dt hths0{t - s) ds 







3=0 



2 ;-T 
f I 

"'0 



where the last equality was obtained by expressing hu = Jq hgds in the second term and changing 
the order of integration over s and u. Expansion of the right member of eq. ( 60 1 gives in turn 



Jexp{-(3TH{x))dx ,^(Rh ^2/rr2\ 



1 + Po ht{H{xt))odt + 



PI 



2 pT 



dt hth,{H{xtf)^ds. 
Jo 



The comparison of the terms quadratic in h leads to the identity 



1 



/3( 



du- 



Js 



3=0 



(61) 



for s < t which gives the relation (59 1 by the derivation with respect to s. Once again the 
Jarzynski equality appears as a global version of the FDT. 



5 Conclusions 



We have discussed fluctuation relations for diffusion processes ([TJ in a non-stationary thermal 
bath. Those included the fluctuation relations for the entropy production (12 1. The work 
performed on the system no longer verifies such fluctuation relations, but that there still exists a 
relation (28 1 which permits to extract the free energy difference in a non-equilibrium experiment. 
We proved that the fluctuation relations involving the functional (34 1 are global versions of the 
Modified Fluctuation-Dissipation Theorem (MFDT) (42 1 around a non-equilibrium diffusion 
extending the MFDT obtained before in [6] and of the usual FDT for energy (|59j |42] 
resulting from a pulse of temperature. On the way, in Sec 4.1 we illustrated the extended MFDT 
on a simple example of a harmonic oscillator in a thermal bath with variable temperature and 
we investigated the physical meaning of the effective temperature introduced in |12l |9] for such 
a system. One should underline that the interaction with a non-stationary bath is one among 
many ways to thermally drive a system. For example the thermodiffusion effect (or Sorret effect) 
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which appears in a bath with non-uniform temperature (i.e (3{x) — has been explained in 

|49] . for the unidimensional case, using a stationary non-equihbrium microscopic model of the 
type ([l]) with G — but with the thermophoretic force — added to the drift (jiji. In the 
same spirit, many years ago, Landauer [30] proposed a model with the wall temperature varying 
along a very narrow pipe filled with the Knudsen gas described by a stationary non-equilibrium 
microscopic model of the type ([l]) but with the thermophoretic force — and the chemical 
force —T^ added to the drift. Finally, other way to drive a system is to consider fluctuating 
coefficients in ([ij, for example [l] considered a fluctuating mass and |44l l33] |48] a stochastic 
friction. It would be interesting to describe the fluctuation relations in those setups. 
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